Bloch oscillations of quasispin polaritons in a magneto-optically controlled atomic 
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We consider the propagation of a quantized polarized light in a magneto-optically manipulated 
atomic ensemble with a tripod configuration. Polariton formalism is applied when the medium 
is subjected to a washboard magnetic field under electromagnetically induced transparency. The 
dark-state polariton with multiple components is achieved. We analyze quantum dynamics of the 
dark-state polariton by some experiment data from rubidium D 1-line. It is found that one compo- 
nent propagates freely, however the wavepacket trajectory of the other component performs Bloch 
oscillations. 
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I. INTRODUCTION 



Quasiparticles are regarded as collective excitations of 
many elementary particles, as well as the mixtures of 
different elementary excitations, which are the basic con- 
structions of matters together with elementary particle. 
Quasiparticles are crucial for understanding many phe- 
nomena in condensed matter physics [l|, Q • One of the 
fundamental phenomena in condensed matter physics is 
the dynamics of particles in a periodic potential under 
the influence of a static force. As it is well-known that a 
quantum particle in a periodic potential possesses en- 
ergy eigenvalues forming Bloch bands and delocalized 
eigenstates known as Bloch functions, the particle under- 
goes uniform motion. After a static force is applied, the 
eigenstates become localized, and the energy spectrum 
becomes discrete with the formation of Wannier-Stark 
ladders 0-|5|. Contrary to common sense, this particle 
oscillates instead of infinitely accelerating by the force, 
i.e., the famous Bloch oscillation d, 0]- In physics, such 
an oscillation is generally explained as a Bragg reflection 
of the accelerated particle, which causes a wave packet 
to oscillate rather than translate through the lattice. 

With the advances of atomic physics and quantum 
optics over the last decades, considerable attention has 
been paid on quasiparticles of light-matter interaction 
since they are suggested as new systems to simulate a 
variety of many-body phenomena 'S, '9] with unprece- 
dented precision and control. A prominent phenomenon 
in light-matter interaction is electromagnetically induced 
transparency (EIT), where the transmission of a probe 
beam through an optical dense medium is manipulated 
by means of a control beam [13, [HI . Stimulated by the 
construction of both quantum memories and quantum 
carriers free of quantum decoherence, dark state polari- 
tons emerge in storing and transferring quantum state of 
light to collective atomic excitations of matter in EIT. 
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A dark state polariton (DSP) is a quasiparticle, which 
is a bosoniclike collective excitation of a photon and an 
atomic spin wave [HI, [l^ . The particle nature of dark 
state polaritons, possessing an effective magnetic mo- 
ment, has been demonstrated by the enhanced deflec- 
tion of the laser beam after light propagates through a 
A-type atomic vapor with a magnetic-field gradient ap- 
plied [11 [l5| . and an effective Schrodinger equation is 
derived to exhibit the wave-particle duality of the dark 
state polaritons [l^ . Bloch oscillations of single DSP are 
proposed [l^l, and a method is described to create an 
effective gauge potential for a DSP [isj . 

Another remarkable property of DSPs in EIT is the 
presence of multiple "spin" components, which open up 
the possibility to study a variety of many-particle ef- 
fects in effective magnetic fields. A phenomenon of bire- 
fringence in EIT has been predicted as a generalized 
Stern-Gerlach effect of quantized polarized light [l9l[20|. 
Collapses and Revivals of DSPs are also experimentally 
demonstrated [2l|, Here, we study spatial motion 

of the DSP with two components in two inhomogeneous 
magnetic fields consisting of a periodical magnetic field 
and a magnetic field gradient. The magnitudes of the 
magnetic fields vary in a direction transverse to the in- 
cident direction of the probe beam. This investigation is 
inspired by the following technical advantages of DSPs: 
refractive index modulation straightforwardly creates a 
scalar potential; a direct measurement is simple for pho- 
tons; and the waveguide and resonator techniques can be 
used to confine the spatial motion of polaritons to lower 
dimensions. In this paper, the DSP with two compo- 
nents is generated in an atomic ensemble with a tripod 
configuration, which controlled by a specially designed 
magneto-optical manipulation based on the EIT mech- 
anism. After obtaining the equation of motion for this 
quasispin DSP by the perturbation theory, we employ the 
single-band and tight-binding model to give an analytic 
treatment of the dynamics of DSP. It is found that one 
DSP component acts as a free particle, the other com- 
ponent experiences a harmonic motion with angular fre- 
quency proportional to the steady force and its momen- 
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We now introduce components Ej (r, of the probe 
field which vary slowly in space and time p5| 



and slowly varying operators 



FIG. 1: (Color online) Tripod scheme of matter- light interac- 
tion involving a linear-polarized probe field and one control 
field in the presence of a magnetic field parallel to the field 
propagation direction. 

turn growing linearly in time (i.e. Bloch oscillations). 
The experimental feasibility of the proposed scheme is 
also discussed using rubidium 87 Dl-hne. 

This paper is organized as follows. Section HIl describes 
the theoretical model for four-level atoms with a tripod 
configuration in external fields. In Sec. IIIIl the quasispin 
DSP has been introduced as dressed fields to describe 
the spatial motion of collective excitation. Afterward, we 
present the evolution equation for this quasispin DSP. In 
Sec. lIVl we present the results of detailed calculations for 
the propagation of the DSP in a washboard magnetic field 
within single-band and tight-binding approximation. We 
conclude our paper in the final section. 



aes = exp(ifccz), (2a) 
o-ej = (7ejexp(ifcz),(j = 1,2). (2b) 

Here continuous atomic-flip ap^i, (r, t) = 
{i/Nr)J2rjeNr'^ii'(^) averaged over a smaU 

but macroscopic volume containing many atoms 
Nr = {2N/V)dV » 1 around position r, where 2N 
is the total number of atoms, V is the volume of the 
medium [2^ [2^, and af^^ (t) = is the flip 

operator of the jth atom. Above and hereafter, we 
take h = 1. Under the electric-dipole approximation 
and the rotating-wave approximation, the interaction 
Hamiltonian reads 

2N f 

Hi = — d^r[{Siaii + S2(J22 + ScCTss) 

- {flcTes + gEiGel + gE2(Te2 + h.C.)] , (3) 



II. THEORETICAL MODEL FOR TRIPOD 
ATOMIC ENSEMBLE IN EXTERNAL FIELDS. 

We deal with an ensemble of 2N identical and noninter- 
acting atoms. Each atom is characterized by three lower 
states {|1), |2), |s)}, and an excited state |e), as schemat- 
ically shown in Fig. [TJ The two ground states |1) and 
1 2) are the degenerate Zeeman sublevels. Atoms inter- 
act with a linear-polarized probe field with frequency v 
and wave vector k propagating in the positive z direc- 
tion, and a classical control light with frequency Vc and 
wave number kc. Here k and kc are the wave numbers to 
the central frequencies v and Vc of the probe and control 
fields, respectively. The excited state |e) with to|, = 
is coupled to the ground state |1) (|2)) with Wp = — 1 
(m|- = 1) via a a+ (ct") component Ei{E2) [HQl- The 
control field with Rabi frequencies Vt is assumed to be 
uniform throughout the sample, and it drives the transi- 
tion \s) — !■ |e) with detuning 5c, which create transparent 
windows for two components of the probe field. We as- 
sume that the detuning of component Ei{E2) from its 
corresponding atomic transition is caused by the applied 
magnetic fields along the z axis with the corresponding 
amount 5i = jitB, the magnetic moments Hi = uipgp^iB 
are determined by the Bohr magneton fi b , the gyromag- 
netic factor gp and the magnetic quantum number rup 
of the corresponding state 



in a frame rotating with respect to the probe and driving 
fields. The parameter g = {e\d\l) y/f / {2eoV) character- 
izes the strength of coupling between the probe field and 
the atoms. Due to the symmetry of the states |1) and |2), 
the transition matrix elements in the above equation are 
equal. Although we assume three-photon resonance in 
the absence of the magnetic field, that is, v — uj'^i — 
Vc — Wgj, where lJ^i — '-^ti ^-nd uP^s ^'^'^ the atomic reso- 
nance frequencies without the applied magnetic field, we 
note that the interaction Hamiltonian Hi is the general 
expression in the interaction picture for a tripodlike link- 
age pattern. Here, (i = 1, 2, c) is the detuning between 
the laser field and its corresponding atomic transition. 



III. EVOLUTION EQUATION FOR THE DSP 
WITH TWO COMPONENTS. 

When a light pulse enters a medium, photons interact 
with atoms of the medium. They are converted into com- 
posite quasiparticles of the radiation and atomic excita- 
tions known as polaritons. In EIT system, there are two 
types of polaritons — - the DSP and the bright state po- 
lariton. To obtain the equation that describes the dark- 
state polariton, first we deal with the dynamics of this 
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atomic ensemble. The Heisenberg equations read 



0'21 
0'e2 



[i ((52 - ^i) - 7c] 0-21 + ig*El (726 - ig*E} o-ei,(4a) 

-iSl - (Tel - ig*E^ {ail - (Jee) " 19* (721 

-«r2*cr^2, (4b) 
[i {6c - 6i) ~ 7c] (Tsl + ig*E^ (Tse - if^Vei, (4c) 

-1^2 - 0'e2 - ig*E^ ((722 ^ CTee) - ig*-^!^ (7l2 

-^^^V,2, (4d) 
[i ((5c - (^2) - 7c] crs2 + ig*E2crse ~ i^^e2, (4e) 

r 



(4f) 



z(5c - - j CT^e + ig*Eiasi - ig*E2as2 

— ifl (fJee — <7ss) , 



where we have phenomenologically introduced the coher- 
ence relaxation rate 7c of the ground states and the decay 
rate T of the excited state. Here, we consider the case 
where the intensity of the quantum field is much weaker 
than that of the coupling field, and the number of pho- 
tons contained in the signal pulse is much less than the 
number of atoms in the sample. Therefore, the perturba- 
tion approach can be applied to the Heisenberg equations 
of the atomic part of the order in gEj , which is introduced 
in terms of perturbation expansion 



= <^^°J + ^^li^ + • • • 



(5) 



In the above equation, /i, G {1, 2, s, e} and A is a con- 
tinuously varying parameter ranging from zero to unity. 



Here, aj^J is of the zeroth order in gEj and altJ is of the 
first order in gEj, and so on. We retain only terms up to 
the first order in the signal field amplitude since the linear 
optical response theory can sufficiently reflect the main 
physical features of the spatial motion of the input pulse 
with slow group velocity. With the assumption that all 
atoms are initially in level {j = 1,2) without polar- 
ization, i.e., atom i is in a mix state pi = I]j=i.2|i)(j|/2, 
we can neglect the population of states |e) and |s), as weU 
as the coherence between these states. Then the first or- 
der atomic transition operator a^j ^ satisfies the following 
equation 



(1) 



(T 



(6) 



which is related to the atomic linear response to the probe 
field. In the above equation, 

dj = « {Sj - Sc) - 7c, 

The dark-state and bright-state polaritons are de- 
scribed by the field operators as 



cos 61 -2%/iVcr,^P sin 6*, 



(7a) 



$j = Ej sine + 2VNal]^ cose, (7b) 



respectively, which are superpositions of photonic and 
spin-wave excitations. They are introduced in the linear 
regime with respect to the probe field. Here, tan 9 = 
g\/N /fl. The Heisenberg equation for the slowly varying 
field operator Ej (r, t) results in a paraxial wave equation 
in classical optics 



■ ^ C , ^2 Q2 ■ 



E,=-2g*Naf^. 



(8) 



Here, c is the velocity of light in vacuum. In terms of 
dark and bright polariton field operators, equations (|6]) 
and read [H 111] 



idt + icd, 
g*VN 



— {dl + d^y)\ (*, cos0 + ci>,sin( 



n* 



{dt - dj) ($j cos 9 - *j sin ( 



(9) 



Under the condition of EIT, absorption is negligible, the 
excitation of the bright-state polariton field $j vanishes 
approximately. Then the dynamics of dark-state polari- 
ton field j is obtained. To show the basic principle of 
physics, we assume a sufficiently strong driving field such 
that |r2|^ ^ r7c. In addition, we have let 7c = 0. To 
make an analogy to the schrodinger equation, we rewrite 
the dynamic equation for the two dark-state polaritons 
as 



p2 



p2 

y 



2m 



+ fif^B{x) 



(10) 



Here, effective magnetic moments 

fJ'f-'^ = (/^s - fij)sin^9, 

which can be adjusted by the control field. In addition, 
the group velocity Vg = c cos^ 9 along the z direction can 
be controlled by the amplitude of the control field, m = 
k/vg is the effective mass, and Pa ~ —ida{a = x,y,z) 
is the momentum operator on the a direction [l9j . The 
steady atomic response in the applied magnetic fields in- 
duces an effective potential for the DSP. Equation (fTO|) 
indicates that the two components of a DSP behave in- 
dependently, i.e. two components of a DSP does not 
interact with each other. Now, we assume that the mag- 
netic fields in z direction consist of two components, one 
is linear with the expression Bix and the other is peri- 
odical in the transverse x direction. Then the treatment 
proceeds along the x — z plane. The evolution equation 
for the quasispin DSP is rewritten as 



VgPz 



p2 

2m 



F,x + V,{x) 



J' 



which means each component j is subject to a static 
force Fj = yf^^^ Bi in its corresponding periodic potential 
Vj{x + d) = Vj{x). Since the effective Hamiltonian along 
z direction commutes with that along the transverse di- 
rection which refers to Wannier-Stark Hamiltonian Q, 
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FIG. 2: (Color online) A possible experimental realization of 
the tripod scheme in matter-light system in rubidium 87 (a) 
and an illustration of the periodical potential for the second 
component of DSP induced by the periodical magnetic field 
(b). 



therefore, the investigation is confined to the transverse 
direction. We note that the evolution operator along z- 
axis is equivalent to the translation operator, and the 
Hermitian of the Wannier- Stark Hamiltonian preserves 
the number of the particle along x direction. 



IV. QUANTUM DYNAMICS OF A QUASISPIN 
DARK STATE POLARITON 

To show the possibility that the dynamics of DSP 
mirrors the Bloch oscillation in this atom-photon sys- 
tem, we consider the atomic medium with light tuned to 
the rubidium (^''Rb) Dl-line 5 S^^^ ^ P1/2 shown 
in Fig. [2Ia). The ground states contain two hyperfine 
ground levels with F = 1 and F = 2, where |1) and |2) 
correspond to the magnetic sublevels (with mp = 1 and 
— 1) of the F — 1 hyperfine ground state and |s) repre- 
sents the hyperfine ground state \F — 2, mp = 1). In this 
case, the effective magnetic moment ^^'^^ vanishes due to 
/is = /ii , which leads to the free propagation of the com- 
ponent ^'i. Therefore, a wave packet with a Gaussian 
profile in space 



1 



(11) 



at the initial time, will keep its shape in all directions. 
Since equation (jlip has a pronounced peak with width ax 
situated at moment px = 0, the wave-packet will always 
localize around the position x — along x-axis, but the 
center of the wave-packet in z-direction moves to z = Vgt. 

Different from 'i'l, the second component ^'2 feels po- 
tentials due to its nonzero effective magnetic moment. 



Obviously, in the absence of the linear magnetic field, 
the period potential gives rise to the Bloch bands. To 
make this discussion more quantitative, we assume that 
V2 {x) is the Kroning-Penney potential which is formed 
by a periodic sequence of rectangular wells with am- 
plitude Vo = 79.15kHz. As illustrated in Fig. Mh), 
the lattice constant of the periodic array d = 8/im is 
equal to a + b. Within the unit cell, the barrier region 
< a; < a, the well region —b<x<0. By tak- 
ing the effective transverse mass m = 7.9 x lO^sm"^ 
for the wavelength of the probe field A = 795nm, we 
obtain the width A = TA.QkHz of the ground energy 
band, and the gap interval Eg = 266kHz between two 
lower bands. For a linear magnetic field with magni- 
tude Bi — 8.5 X lO^/iGmm^^, the magnetic moments 
= -4.64 X 10-24 jy-i and /i, = 4.64 x 10^24 jy-i 
allow us to find that F2d Eg. So the ground band is 
well separated and one can neglect the interband mixing 
induced by the static force F2. By assuming that neigh- 
boring wells are directly coupled, the dynamics of the 
DSP component ^'2 is described by the Hamiltonian 

Hi = ^^(|,i) {n + l\ + h.c.) + dF2j2^\^) 

n n 

(12) 

in terms of Wannier state which is localized around 
position X — nd. It is useful to introduce the eigenstates 
\k) (Bloch states) of Hamiltonian Hi with F2 = 0, 



2n ^ 



(13) 



which is the Fourier transform of the Wannier states. The 
tight-binding Hamiltonian in Ea. (|12p is diagonalized as 



(fc) = cos {nd) -t- iF2^ 
2 dn 



(14) 



in the quasi-momentum representation with the quasi- 
momentum k, G [— Tr/d, Tr/d] in the first Brillouin zone. 
Diagonalizing Hi gives the Wannier-Stark ladder with 
energies Em ~ mdF2, where m is an integer. The 
Wannier-Stark state \ipm) to the eigenvalue Em is given 

by 



\4'n 




^^g-i[A sin{Kd)/{2dF2)+mdK] 



In the Wannier-state representation, we obtain 



^ ^ Jn — rn f 



\2dF2 



In) 



(15) 



(16) 



where Jn-m is the Bessel functions of the first kind and 
of integer order n — m. 

To show the wave-particle duality of the second DSP 
component, we now discuss the quantum dynamics under 
the assumption that an initially Gaussian wave packet 



l'^(O)) 



1/4 



exp 



n^d^ 
1^ 



(17) 
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is launched into the lattice with a center value at Xq = 
and width a. Since the time-evolution operator is di- 
agonal in Wannier-Stark states, each element [/^^ = 
U (t) Iv) of the time-evolution operator U (t) in the 
Wannier-state representation has the form 

U^, = J,_^ (^1^) e-'^-..*e-'^(''-'^)K*--)/2, (18) 

where a = ujgt/2 and — dF^. The argument of Ji/-^ 
oscillates in time. The state at arbitrary time can be 
expressed as — /« (^) l*^)- Gaussian 

wave packet is broad initially, i.e. a /d^ 1, the coefficient 
/„ [t) is approximately given by 

where rit — Xc (t) / d. Here, the motion of the wave packet 
center along x-direction 

a^c(i) = ^[cosKi)-l] (20) 

reduces the quantum mechanical dynamics to the clas- 
sical trajectories. Equation ([20)1 implies that the DSP 
component ^'2 experiences a harmonic motion around 
the initial position xo with angular frequency lo^ and 
amplitude A = IS./ {2F2). As expected from the semi- 
classical picture, the band width A equals the product 
of the total spatial extension 2 A and the static force F2. 
Besides the trajectory of the wave packet, equation (fT9|) 
also gives a description on the variation of the particle's 
wavenumber k = uj^t/d ~ F^t in the semiclassical pic- 
ture. Since z = Vgt is the relation between the center of 
the wave packet along the z direction and the time, we 
achieve the spatial Bloch oscillation 

a; (z) = ^ [1 - cos {(^z)] 

with spatial period 2Tr/(. Here, C = dF^/vg. In Fig. [3l 
we plot the probabilities for the components Vl/i and 4*2 
versus x/d and i, where the initial states of both com- 
ponents are identical broad Gaussian wave packets. One 
can see that two components behave differently in the 
transverse direction though the wave packets for both 
components keeps it shape along z direction with an un- 
changed variance. The trajectory of component ^1 is 
always localized at its initial position, which means this 
component appears like a free particle. However, com- 
ponent ^'2 presents a back and forth behavior around an 
equilibrium point with a constant amplitude and a con- 
stant frequency, so it acts as a particle moving through 
a quadratic potential. 

V. CONCLUSION 

We study an optical-controlled cloud of identical atoms 
in a washboard magnetic field applied along z-direction. 
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FIG. 3: (Color online) Density plots of the norm square of the 
wave function for a DSP with two components arising from 
an extended Gaussian distribution. The probabilities for the 
components ^'2 (a) and 'I'l (b) is plot as a function of x in 
units of lattice constant d and t in units of millisecond, where 
the initial states for both components are taken to be the 
same with width 0.1mm. 



The probe and control laser beams drive three hyperfine 
ground states to a common excited state, which form a 
tripod configuration. The quasispin DSP is excited in 
the EIT condition. The equation for the space-time evo- 
lution of this quasispin DSP shows that each component 
is subject to its corresponding effective potential induced 
by the steady atomic response in the external spatial- 
dependent field. By taking the value of parameters from 
the experiment data in rubidium D 1-line, it is found that 
one component of the DSP acts as a free particle. The 
behavior of the other component '^2 is analyzed via the 
single band and tight-binding approximation. By con- 
sidering the time evolution of the broad Gaussian state 
for the component 4*2, the particle nature of component 
^'2 is shown by the trajectory of the wavepacket, which 
undergoes a periodic motion with angular frequency pro- 
portional to the steady force. Besides, its momentum 
grows linearly in time. This periodic motion is termed as 
Bloch oscillations. The oscillation amplitude and period 
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are controlled by the intensity of the control field and the 
magnitude of the magnetic fields. 

It should be pointed out that there are many differ- 
ences between our present scheme and that in Rcf. 
1) The atomic configurations are different. Lambda type 
is studied in Ref. [17|, however tripod pattern is studied 
in our manuscript. From the technical point of view, tri- 
pod atoms proved to be robust systems for "engineering" 
arbitrary coherent superpositions of atomic states [27[ 
using an extension of the well known technique of stimu- 
lated Raman adiabatic passage. 2) The double-EIT effect 
has been used in our manuscript. Double EIT is an im- 
portant phenomenon for quantum information processing 
and quantum computing. 3) In our paper, different types 
of dark polaritons are demonstrated to possess different 
effective magnetic moments. 4) The predicted phenom- 
ena in our manuscript are more general. It could be 
found that only one dark polariton is excited in Ref. [13] , 
so only one trajectory is found inside the atomic medium. 
However, in our paper, two dark polaritons are excited 
by one probe beam and a spatial bifurcation and dynam- 
ics of these dark polaritons is obtained inside the atomic 
medium. 5) Our approach naturally shows the multiple 
degrees of freedom of photons and the role of quantum 



coherence. As one can find that when the single pho- 
ton with a superposition of two orthogonal components 
of polarization is incident on the atomic medium, one 
obtain a superposition state of two DSPs relating to dif- 
ferent spin-states of this incident single photon. Besides, 
our investigation provides a way to control the dynam- 
ics of DSPs. In addition to the external fields, atomic 
energy levels can be used as a way to adjust the motion 
of DSP, for example, one can choose the atomic energy 
level to make the DSP have vanished effective magnetic 
moment. The presence of the first component unaffected 
by the washboard potential might be used to eliminate 
the distortions or aberrations in future. 
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